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On the determination of the laser modulation index from the optical 
absorption spectrum for clock implementation by CPT 



Introduction 



In the use of the CPT transmission absorption line for implementing a frequency standard it is 
essential that the laser modulation index be set at a given value. In particular, it is interesting • o 
set it at m = 2.4, the value that makes the power light shift equal to zero. In that case the clock 
frequency 1S essentially mdependent of laser power. It can also be set at about 2.6-2 8 where 
diere is a minimum » .the .light shift curve against modulation index, making the clock frequency 
independent of the index in first order. Which setting is best regarding frequency stability is not 
yet known. Experiments will decide. «Miuy»noi 

However, before adopting either approach, it is necessary to see if we can set easily the 
modulation indexto a prescribed value without too much effort. In the standard method a 
Fabry-Perot interferometer is used for that purpose. It is the first tool that comes to mind 
However the setting of the instrument and the measurement process is complex. Furthermore it 
requires the laser to be examined independently of the actual clock setup. An approach in which 
the laser modulation can be evaluated in situ without disassembling the clock should provide a 
net gain. Such an approach is possible through the direct examination of the light absorption 
spectrum at the exit of the cell by means of the same photodetector diode as that used in the 
drtection system of the CPT transmission resonance line. This is what is examined in the present 

I^SfT °^ he modu ! ated J aser is calculated at the exit of the Rb gas cell, as observed at the 
photodetector. The analysis is done as a function of the modulation index of the laser The 
analysis is done also as a function of light intensity, Rb density and buffer gas pressure The 
results are compared to the experimental data obtained in the laboratory. Conclusions are drawn 
on the feasibility of the technique. 

The analysis does not take into account optical pumping taking place at high light intensity 
However, smce the theoretical results neglecting this effect are in good agreement with the 
expenmental data, it is concluded that the technique can be used in practice with reasonable 
accuracy. 



A future report will discuss in more detail the effect of optical pumping and provide additional 
modSo°n 00 ^ Ue " Ce m the im P lement ation of the technique for adjusting the index of 
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It should be realized that the technique could be included directly in the start sequence of a 
frequency standard using the digital approach. 



The setup 

The standard CPT setup is shown in Fig .1 . 
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Figure 1. Experimental setup. 

The output of the transimpedance amplifier is V d proportional to the light intensity falling on the 
detector. V d is a measure of the total light spectrum. There is no frequency discrimination in the 
detection process and the output voltage is the sum of all the laser sidebands after absorption by 



The absorption spectrum. 



The Rb atoms in the cell are characterized by the energy level structure shown in Fig. 2. 
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Figure 2. Illustration of the four hyperfine transitions of the Dl line of 87 Rb with their respective 

transition probability. 

If the laser is monochromatic, not modulated, the cell acts as a frequency discriminator providing 
four absorption hues a, b, c and d, easily identified by sweeping the laser frequency slowly over 
foe spectrum. The W1 dth of these lines is of the order of 700 to 1000 MHz depending on the 
buffer gas pressure. The hyperfine lines are all resolved, although overlapping to a certain 
extent. A typical spectrum is shown in Fig. 3 with the respective lines identified. It is readily 
observed that line c has a much lower intensity than the others do, as predicted by the transition 
probability shown in figure 2. 
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Figure 3 Typical absorption spectrum observed in 87 Rb in buffer gas. 

It is observed that when circular polarization cr + or <T is used, there is reduction of the 
absorption of line a The most likely explanation is optical pumping to a Zeeman level that is not 
in interaction with the laser radiation. The population of the absorbing level is reduced and the 
absorption involving that level is also reduced. When the laser frequency is scanned across the 
absorption spectrum there is also optical pumping from one hyperfine level to the other hyperfine 
level. If the scan is slow as is normally the case, an equilibrium is reached where one ground 
state hyperfine level is more populated than the other, reducing the optical absorption. The first 
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process is present in the observation of the dark line because circularly polarized light is required 
for the observation of the 0-0 transition. The second process is not present to the same extent in 
the observation of the dark line since the two ground state hyperfine levels are equally pumped 
by the two radiation fields J I+ and J,, of equal intensities. Experiments show, as expected that 
optical pumping is more important at high light intensity. 

These effects will be examined in more details in a future report and their consequence on the 
present analysis will be outlined. 

Theoretical background 

The radiation amplitude of the "n"th sideband in the laser spectrum is described by the electric 
field Eo,,. We define the Rabi frequency proportional to this electric field as : 

©Rnij = (Eon/rT) <i |er«e x |j> ( 1 ) 

This definition is introduced in order to simplify notation and provide better insight into the 
physical mechanisms taking place in the laser radiation absorption process. In that equation, n is 
the sideband identification, tT is Planck's constant over 2*. and the terms between brackets 
represent the electric dipole matrix element characterizing the transition between levels i and j 
It is generally written as dy and gives the intensity of absorption. 

Absorption is described by the differential equation derived from the Maxwell's field equation 
coupling the radiation electric field to the polarization of the Rb ensemble. The polarization of 
the Rb ensemble is calculated in the density matrix formalism through solving the appropriate 
rate equations for the level populations and the coherence existing in the system and introduced 
by the laser radiation. For sideband n and transitions between levels i and j an exact calculation 
gives: 

-^- =a « lm5 nij (2) 

where a is the absorption coefficient defined as 



«y = 



v ce 0 /i "j 



n Rb (3) 



All the effects of optical pumping and coherent population trapping are imbedded into the term 
ImSnij, which means imaginary part of the off diagonal density matrix element o njj . It is the 
optical coherence created in the system by the radiation field sideband E„ at the transition 
frequency coresponding to the transition between levels i and j . The transition probability for 
transition i to j is imbedded in the matrix dipole moment dy. On the other hand, the various 
terms in Oy are defined as follows: a> is the average laser frequency, c is the speed of light, e„ is 
the permittivity of free space and n Rb if the Rb density. 

If we neglect optical pumping from one level to another level of the ground state, ImS,,;; is given 
by; 



Im 8 nij = - 



Ur/2) 2 + <fi nlj) *J 



(4) 
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where Q ni j is 



f2nij - CD„ - G)jj 



(5) 



ta!l&L te ,aser sideband angular frequency and ^ angular frequenc y of * c atomic 

In the theory parameter r is the decay rate from the excited state caused by Rb-buffer eas atom 
colhaons. Unfortunately, there is always broadening from Doppler effect Ld in pSti^C 

2h T ^ ? ^ ^ ** CXpeCted JUSt fr0m the e * cited decayme Actally 
Uie optical absorption line ,s a convolution of a Gaussian line shape (Doppler effect) and of a 

„TLt^ ( ??• CXCited StatC: V ° igt profi,e )- In that the problem is 

intractable since the solution of the above differential equation would need to be integrated Tover 

all velocit.es. However, smce in practice the line shape observed is closely LorenteiaTttis 
possible to approximate die situation by assuming a decay rate that gives an absorption line 
width the same as the one observed. This is the approach we use. In that case the differential 
equation can be integrated directly and gives Beer's law for absorption- 



(0 



Rn 



< z ) = °>Rn(0)exp -a 



(174) 



Ur/2) 2 + (« (jn ) : 



(6) 



^iS 6 "^ fate 8iVing 3 Widdl CqUaI t0 (1/2n)P ' ^--ating the 
In this expression 0^(0) is the value of the Rabi frequency at the entrance of the cell. According 
to Eq. 1 it is proportional to the radiation electric field of the nth sideband. The voltage 
measured at the photodetector of the setup shown in Fig. 1 is proportional to the intensity of the 
radiation, thus to the square of the electric field of the radiation. Furthermore this voltoge is 
proportional to the sum of all the radiation fields traversing the absorption cell, that is all the 
sidebands. Consequently a summation must be made over all these sidebands n. Furthermore a 
summation must also be made as well on all the absorption lines <i[j> shown in Fig 2. The result 

IS - 



K(z)) 2 = I (o)^ (O)) 2 exp -21 



r 



074) 



| v (r/2) 2 + (o IJn ) ; 



) 



(7) 



We have also introduced the coefficient a* that takes into account the actual transition probability 
shown in Fig. 2 and leaves a as a general term constant for all transitions. 

Since V d is proportional to the square of the Rabi frequency this equation can be written as 

V d = kl < ttRn (0)) 2 exp .21 A <£/£) 

« 'l( r/2 ) 2 +^ijn) 2 J 
Here k is a constant representing the transformation of light intensity (Rabi frequency) into 
voltage by the detection system. 



(8) 
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Approximations made 



The constant to be used 



s^htT^e £iM^^ l S P ^r B diSCUSS6d ab0 - * P-ticc it is set 

line shape ITie^ueC^l fo^a eel , T eX P erimenta "y. turning a Lorentz 

10V. value used here for a cell containing a N 2 -Ar buffer gas mixture at 10Torris4x 

«tssSMSs: 2?cS3r caicu,ation - *• contrast " S sion 

expeiaidatawi^l^ 

« ■ » " takCn 38 ^ *™ * *■ « ^ 1 for a„ transitions and 0, for 

The value oj 'the Rabi frequency at the entrance of the cell <o Rn (0) We set it for thm • r 
unmodulated laser. We assume a value eaual to 2 * 1 o 6 t2 I r ! he Camer > for 3,1 

then obtained through a multi D licatim h^A-a • I * f ° r the Various side b*nds is 
modulation chosen mulU P hcatlon b y *e appropriate Bessel function value for the index of 



The calculation. 



£2^"^^ ^ «■* Ch ° Sen ab ° Ve - A- -ample of 
Fig. 4. Detailed results^ f «,ual to 1.8 is sho™ in 

calculation. y 1,16 Jz> J| 31,(1 J ° sidebands are used in the 

ra=1.8 w R =Intensity of radiation transmitted W=Afrequency 
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Figure 4 Theoretical absorption spectrum of 87Rb, Dl resonance optical line. The laser is 
assumed monochromahc and modulated with an index of 1.8. The constants used Z 'the 'model 

are explained in the text. 



Experimental results 



A typical experimental result for a modulation index of -1.5 is shown in Figure 5. 



6 O H- 79 6 8 7 . O 6. JL 9 O 



Agilent Technologies 

a sow 




* 




Figure 5. Experimental absorption spectrum in " Rb observed on an operating CPT Clock with 

frequency modulated VCSEL. 



to evaluate the index of modulation to approximatefyTo l -ST » sufficient 



Determination of the index of modulation 



The index of modulation can readily be evaluated by plotting the ratios flusn and rci /<w \ 
These terms are defined in Fig 4. Due to the form ofEq. 7, fhese raUo 2e « enSuy r fios of 
Bessel fonctmns included in the values of These ratios are plotted for the Scd resuL 
m F,g. 6. Expenmental results for the three optica, densities studied are ^SSS^Tt 



Theoretical data mod Index 
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Figure 6. Theoretical ratio of resonance and satellite lines 



Ratio of sat. abs. Lines OD 0.3 



♦ Rx/Six aSlx/S2x 




Power [mW] 



200 



250 



Figure 7 ^mental data on resonance and satellite ratios for light intensity corresponding to a 
gray filter with optical density equal to 0.3 placed in the path of the laser beam 
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Figure 8. Experimental data on satellite ratios for light intensity corresponding to a gray filter 
with opucal density equal to 0.5 placed in the path of the ikser ocam 
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ntodulation «„ better than IM by Zp * rn^Th. S * or S ?S **" ? f 
125 mW, 21 dBm, the index is about 2.4. " " S2 ' For example 31 

At higher power ( >150mW), it appears that there is saturation in the system It is not clear 
where tjtts effee, contea fl™. „ is poss i W e ttta, the synthesizer is JETJ£E5L 
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Effect of optical pumping rate, line width and temperature. 

The same calculation were done for three other situations in relation to the absorption 

increased. This ,s expected since the light intensity does not alter the opdea teS £d 
optical pumping has been neglected. It is also seen that values of r aj i „ „n,„?i » . 
above in the calculation do no, provide better ag^'bt^exSm"*^ 



Conclusion 



It appears that the theoretical results are in fairly good agreement with the experimental data. 

1) Except for small differences in the low frequency absorption lines of the spectra inaht h™* 

ShlK" spectra \r, shapes m rather ■ taB - «- 5* with S3^?3S-2K? 

much like the expenmental data vary with microwave power. The small differences are believed 
to be due to Zeeman optical pumping taking place under circular polarization 

2) Differences between experimental data and theoretical data at high modulation power are 



eX P erimental s P ec <" ■ efficient to determine that 
the experimental case of 22 dBm corresponds to a modulation index m equal to 2 6 while 21 

dBm corresponds rather closely to m = 2.4. Calculation of actual ratio of satellite anSdes 

provides an effecfye mean of determining the modulation index better th^ToSSSS bv 

visual inspection. The ambiguity that results from the quasi-parabolic response of meTe rlnoV 

removed by exammation of the minimum or maximum of the ratios observed. 



Annex I 



sideSsI ' T^nH t ' , T S: f bS ° iptl0n SpeCtmm 0f modulated ,aser radiati °n- Only the 
sidebands Jo, J, and J 2 are taken mto account. J 3 , if considered, would add a small component at 

- , , ,5fft V f UCS ° f modu,ati ? n j nde *- The constant assumed are r = 4 x loV 
a-2.1 xlO m.s and ca R = 2 x loV. These corresponds to a typical situation: T- = 65°C 
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Pbg =10 Torr, and a light intensity that produces and added dark line broadening of the order of 

250 Hz. 

Notation: 

- Muv means an Index of modulation equal to (u.v). For example, raO means that only the carrier 
is present while ml6 means an index of modulation equal to 1.6. 

-The second figure (ml2) identifies lines to be used in the analysis. R t corresponds to the line 
used for dark line resonance observation, while S 1, and S2, are satellites used in the determination 
of the index of modulation. Subscript t stands for theoretical x later used means experimental 
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m-1.5 w R -Intensity of radiation transmitted W=Afrequency 
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■n-1.6 w R =Intensity of radiation transmitted W=Afrequency 
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m1 * 7 ^ R=Intensit y of radiation transmitted W=Afrequency 



4' 10 



3.5 '10 



3*10 



12 



12 



12 



2.5 10 12 



2* 10 



13 



1.5 '10 12 



1*10 



12 



5' 10 



u 




• 



■1 



■ 



200 



400 



-» ■ ~i 1 i_ 

600 



w 



600 



1000 



m=1.8 w R =Intensity of radiation transmitted 
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Figure A2. Absorption spectrum for a modulation index of 2.4 but for three different situations 

regarding absorption coefficient, line width, and Rabi frequency 

Case a: stronger light intensity, a = 2. 1 x 10 1 \r = 4 x 10 9 , a> R =4 x 10 6 , m = 2.4 
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Case b: narrow lines, larger abs. coef, a = 4.5 x 10 n , r = 2 x 10 9 , o R =2x 10 6 , m = 2.4 
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Case c: less absorption but broader line, a = 1.5 x 10" J" = 6 x 10 9 , © R =2 x 10 6 , m = 2.4 
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IMPLEMENTATION SUPPLEMENT TO- 
"On the determination of the laser modulation index from the optical absorption 

spectrum for clock implementation by CPV PATENT 

The functioning of the coherent population trapping frequency standard, as described in 
(reference to Vanier patent number) is based on the interaction of a frequency -modulated 
optical earner with the hyperfine resonance of an alkali-metal vapor contained within a glass- 
windowed cell A typical arrangement is shown in Figure 1 below. The optical carrier is 
generated by a laser diode which is frequency modulated at one-half the hyperfine resonance 
frequency of the alkali-metal vapor. For example, for the Dl line of rubidium87 the hyperfine 
frequency is 6.8 GHz and therefore the laser is modulated at 3 4 GHz nypertine 



Attenuator 



Resonance 

/»nll» 

\ 




Current 

source 

Microwave 
generator 
v = (l/2)v„ r 



A/4 plate 



Solenoid i 



0 ! 



Photodetector 



Signal 




Magnetic shield 




Figure 1 Typical CPT Resonance Configuration 
The frequency modulation process generates a set of predictable sidebands disposed 

Iw?!??* T r/ e,ength - ^ magnitude of these sideba »ds is shown in 

v. a > * ^ ^ modulation index > whi <* is a measure of the degree of modulation 

apphed to the laser source. It is shown in the references that the performance of a CPT-based 
frequency standard is critically dependent on the selection and maintenance of specific values of 
the modulation index. Unfortunately, the experimental determination of the frequency 
modulation index ,s difficult and normally requires a specialized optical spectrum analyzer 
However in the special case of the CPT frequency standard schematically shown in Figure 1 it is 
possible to extract die modulation index from a set of simple measurements and a matoat Lai 
algorithm described in the attached reference by Vanier. mamemaucai 
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Figure 2: Frequency Modulation Sidebands as a Function of Modulation Index 

If the wavelength of an unmodulated laser is slowly sweDt across the hvn«-fi„» 
S, J? resonance features result from transitions Horn the two hyperfine states of 
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Figure 3: Hyperfine Spectrum of Rubidium87 Using an Unmodulated Laser 

However, if the laser source is modulated at approximately one-half the hyperfine separation 
shown in Figure 3 the resulting spectrum traced out by the photodetector current as the laser 
wavelength, together with all of the frequency modulation sidebands, is swept over the hyperfine 
resonances. The spectrum of Figure 4 is mathematically the convolution of the modulated laser 
spectrum with the hyperfine absorption spectrum: a typical pattern is shown in Figure 4 below 
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Figure 4: Hyperfine Spectrum of Rubidium87 Using a Frequency Modulated Laser 

It is shown in the attached paper by Vanier that Figure 4 expresses all the information required to 
extract the modulation index of the modulated laser source. The specific procedure canbe 
summarized as follows: 

a. The laser source is modulated at approximately one-half the hyperfine separation 

b. The center wavelength of the laser is then slowly swept, by mechanical, thermal or 
electrical means, through the hyperfine absorption spectrum of the alkali metal vapor 

c. The output current of a photodetector monitoring the optical transmission through the 
vapor cell is converted to a voltage and digitized by an analog-to-digital converter which 
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is a component of the frequency standard control servo. 

d. The digital values corresponding to each of the minima shown in Figure 4 are passed to a 
microprocessor, which is a component of the frequency standard control servo. 

paper by Vanier, the modulation index 
determined. The modulation power may then be increased or decreased as required 
produce the desired modulation index. 

f. The process is terminated and the functionality of the frequency standard returned to 
normal operation, as described in the references. 
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